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Abstract 

In this paper we have derived the behavior of deceleration parameter with respect to redshift in 
context of f(R) gravity in vacuum using Taylor expansion of derivative of action. Here we have 
obtained that the two first terms in Taylor expansion may describe the late time acceleration which 
is appeared by SNela without need of dark energy and dark matter. Also we have derived that any 
other terms higher than z in Taylor expansion may describe main inflationary epoch in the early 
Universe. We have shown that f(R) gravity may cover all the dynamical history of the Universe 
from the beginning to the late time accelerating phase transition. 
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The recent data coming from the luminosity distance of SuperNovae la (SNela)jT], wide 
galaxy surveys [2] and the anisotropy of cosmic microwave background radiation |3J H] suggest 
that the Universe is undergoing an accelerating expansion. Several candidates, responsible 
for this expansion, have been proposed in the literature, in particular, dark energy models 
and modified gravity. The main problem of dark energy models are understanding their na- 
ture, since they are introduced as ad hoc gravity sources in a well define models of gravity. 
In this context, modified theories of gravity, such as f(R) gravity [3 El [7], has been verified 
in an attempt to explain the late-time accelerated expansion of the Universe. These theories 
are also referred to as 'extended theories of gravity', since they naturally generalize General 
Relativity. There has been predicted that the universe might have appeared from an infla- 
tionary phase in the past. It is also believed that the present universe is passing through 
a phase of the cosmic acceleration. While a series of cosmological models were proposed in 
Einsteins gravity with early inflationary scenario are working well, the mystery of passing 
universe through an accelerated phase of expansion is an interesting open question that is 
yet to be understood. There are many interesting works have done in the context of f(R) 
gravity on unifying late time accelerating behavior and early inflationary behavior of the 
Universe [8], where they have shown that the late time acceleration is in continues of early 
time inflation for specific models of f(R) gravity. Also, their vacuum solutions are obtained 
for constant Ricci scalar while it is possible to derive non constant curvature scalar vacuum 
solutions in f(R) gravity theories. In this letter we show easily that f(R) gravity theories 
are the most powerful approach to describe all the dynamical behavior of the Universe in 
one scenario. Here we look to the story of f(R) gravity from a different angle. Actually we 
consider vacuum solutions of f(R) gravity. But there are two main difference with other 
vacuum solutions. The first one is that we do not assume constant scalar curvature to 
obtain vacuum solution. The second one is that we do not consider any specific model of 
f(R) gravity action. Vacuum solutions of modified f(R) gravity would like to explain all 
the phase transition of cosmological parameters like deceleration parameter without need of 
dark companion of the Universe, alone by pure geometry. The action of modified theory of 
gravity is given by 



where L m is the matter action such as radiation, baryonic matter, dark matter and so on 
which we do not consider in field equation. In this work, we consider the flat Friedmann 
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(1) 



Robertson Walker, (FRW) background, so that the gravitational field equations for f(R) 
modified gravity are provided by the following form 

-3-f' + 3-Rf"+lf = 0, (2) 
a a 2 

[- + 2- 2 ]f - 2-Rf" - R 2 f"> ~ Rf" -b = 0. (3) 
a a 2 a 2 

where the overdot denotes a derivative with respect t, a(t) is the scale factor and H = 
a(t)/a(t) is the Hubble parameter. Eliminating / between Eqs. (j2l) and ^ obtains: 

- 2[- - {-f\f + -Rf" - Rf" - R 2 f" = 0. (4) 
a a a 

Which can be changed to the form of: 

F-HF + 2HF = 0, (5) 

where F = df/dR. Eq. Q is a second order differential equation of F with respect to 
time, in which both of F and H are undefined. The usual method to solve Eq. ^ is based 
on definition of f(R)- But in this paper we would like to replace the variable of Eq. (J5]) 
by another cosmic parameter which named redshift, z. Each redshift, z has an associated 
cosmic time t (the time when objects observed with redshift z emitted their light), so we 
can replace all the differentials with respect to t by z via: 

d dadz d 
dt dt da dz ' 

= -{l + z)H{z)± (6) 

where we use 1+z = ao/a, and we consider ao = 1, in the present time. Now, we can replace 
the variable of Eq. (|5| from t to z by using Eq. ([6]), and we obtain a first order differential 
equation for H 2 with respect to z as: 

±-H{zf = P(z)H(z) 2 , (7) 

where P(z) is a function with respect to z, which depend on definition of F with respect to 
z as: 

2(1 + z)(d 2 F/dz 2 )+4(dF/dz) 
{ ' 2F-(l + z)(dF/dz) ' 1 J 

The deceleration parameter, q in cosmology is a dimensionless measure of the cosmic accel- 
eration of the expansion of space in a FRW universe. It is defined by: 



here we change the variable of Eq. ^ from t to z to have evolution of deceleration parameter 
with respect to redshift as: 

1 + z dH(z) 2 . . 

= -i + ^?W, (ii) 

which only depends on P(z) for vacuum solution. Since we would like to consider a general 
behavior of q(z) due to the model independent F(z), we use Taylor expansion of F(z) around 
the present value of redshift, z = as: 

m 

F{z) = Y,F n z n (12) 



n=0 



where 



1 d n F 

F n 



(13) 

2=0 



n\ dz n 

are some constants at z — 0. Now to calculate deceleration parameter, with some algebra 



we put Taylor expansion of Eq. (12) in Eq. rt8ft, then we put the result in Eq. (11) to find: 



(l + z)^2 nF n z n - 2 [(n + l)z + n-l] 

?(*) = "I + ^ , (14) 

Y.FnZ^^-^z-n] 

n=0 

which is the model independent deceleration parameter in a redshift based f(R) gravity. 
The present time deceleration parameter for a general F(z) obtains as: 

go = q(z = 0) = -1 + (15) 

1 — a/ A 

where a = Fi/F , (3 = F 2 /F which shows that the present value of deceleration parameter 
depends on, only the first two terms coefficients of Taylor expansion of F(z). Since the 
present value of deceleration parameter should be negative, go < 0, a and (3 will constrained 
by 3a/2 < 1 — 0. In the early time, z — > oo the only remained terms in the numerator and 



denominator of the right hand side of Eq. ( 14 ) are the terms with higher powers in z as: 



/ x / s m(m + 1) . . 

qoo {m) = q { z ^oo) = -\ + -± S 16 

2 — m 

which is independent of coefficients of Taylor expansion of F(z). Here we see that 
depends on the number of the sentences in Taylor expansion. For example if m — (one 




FIG. 1. Evolution of deceleration parameter in z — > oo with respect to the number of sentences in 
the Taylor expansion of F(z). 

term in Taylor expansion), we will have a constant modification, which will recover de' Sitter 
space, goo(0) = — 1. If m — 1 (two terms in Taylor expansion), the value of deceleration 
parameter is q'oo(l) = +1, which is consistent with decelerating Universe after inflationary 
epoch. If m = 2 (three terms in Taylor expansion), we have goo (2) — > ±oo, which is an 
unstable value for deceleration parameter. Eventually for any values of m ^ 3, the early 
time value of deceleration parameter is a negative number, q^m ^ 3) < — 1, which shows 
a continues inflationary behavior. One can show that adding any extra term higher than 
first order term in Taylor expansion can make an unstable epoch in evolution of deceleration 
parameter. But location of this instability which we may call it, inflation, depends on the 
value of coefficients of Taylor expansion. Dependency of on the number of sentences in 
Taylor expansion of F(z) is plotted in Fig. ([!]). As a very simple but phenomenological case 
we can assume there are four nonzero terms in Taylor expansion as: 

F(z) = F (l + az + fiz 2 + -fz 3 ), (17) 

where 7 = F^/Fq. Depend on choosing expansion coefficients we can see evolutionary 
behavior of deceleration parameter with respect to redshift. There has plotted three different 
cases for early time period in Figs. (|2||4]). The common result of three above figures is that, 
There is a long era between early time and late time evolution of deceleration parameter 
in which q(z) is equal to +1 in each figure. This period is called radiation dominated era. 
Before that in higher z's there are different instabilities. The other common point is that 
they desire to have an acceleration phase transition point in late time Universe. In Fig. ^ 
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FIG. 2. Evolution of deceleration parameter which is introduced in Eq. (17) for 7 = 0, (3 > solid 
line and (3 < dashed line. 



(3 = ±10 -12 and 7 = 0. For the positive sign (solid line) of /3 the Universe may begin from 
an infinite decelerating state which is continued to the end of radiation dominated era. For 
the negative case (dashed line) the Universe may begin from an infinite accelerating state 
which is continued to the beginning of radiation dominated era. There is one acceleration 
transition point in early time in this latter case. In Fig. ^ we choose /3 — and 7 = ±10~ 25 , 
then for the positive sign (solid line) the Universe may begin from a finite accelerating state. 
But it is faced to an instability in a identified redshift which is called inflation. Therefor 
deceleration parameter goes to —00 and then comes back to a finite value from +00. For the 
negative case (dashed line) the Universe again may begin from the same infinite accelerating 
state which goes to the transition point in a continuous state up to join to beginning of the 
radiation era. In Fig Q there are two cases of composition of values of coefficients. In the 
first one (3 = +10 -12 and 7 = — 10~ 25 (solid line) which shows that the Universe may begin 
from a finite accelerating state and then goes to a finite maximum deceleration value and 
then comes back to radiation era. There could not be an inflationary period in this case. 
While in the other case (3 = — 1CT 12 and 7 = — 10~ 25 (dashed line) there could be an stable 
inflation. Here we saw that smallness of the coefficients f3 and 7 may move any inflationary 
behavior to early Universe. However this example of composition of expansion of coefficients 
are nor realistic but will make this sense that we can find a close combination which might 
be compatible with number of e-folding and other early time cosmological constraints. The 
other main late time cosmological constraint which is tested for the above simple example 



-10 
-15 



10 10 4 10 7 10 10 10 13 

z 



FIG. 3. Evolution of deceleration parameter which is introduced in Eq. (17) for f3 = 0, 7 > solid 
line and 7 < dashed line. 
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FIG. 4. Evolution of deceleration parameter which is introduced in Eq. (17) for the case (5 > 
and 7 < solid line and the case f3 < and 7 < dashed line. 



is SNela distance module with respect to redshift constraint. In this comparison we used 
the Union2 data set [9] which is provide 557 SNela specifications. According to the 
value of numerical results which was taken in comparison with SNela data set from usual 
X 2 algorithm, we obtain that for a = 0.41, h = 0.7 and no need to dark companions of 
the Universe, \ 2 = 546.48 and the ratio of chi square error to the number of freedom is 
about 0.98. The result of this calculation is shown in Fig. ([5]). According to the obtained 
value for of parameter a, evolution of deceleration parameter for low redshift Universe is 
plotted in Fig. Other cosmological parameters such as jerk, snap and lerk is calculated 
for this combination of coefficients. The obtained values of present cosmological parameters 
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FIG. 5. Accordance between Union2 SNela data and late time approximation of Eq. (17) for 
a = 0.41 and h = 0.7. 
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FIG. 6. Late time evolution of deceleration parameter for a = 0.41 and h = 0.7. 

(Qo-ijoi s 0) ^o) are (—0.48, 0.37, —0.03, 0.71) while their values at the beginning of radiation era 
(q r ,j r , s r , l r ) obtain as (1, 3, —15, 105). The same values of present cosmological parameters 
from ACDM for Q M = 0.3 and h = 0.72 is (-0.55, 1.0, -0.35, 3.11). In this letter we have 
considered a general effect of changing the gravitational action on cosmological deceleration 
parameter, from R to f(R)- Here we have obtained that all the cosmic history may describe 
by knowing effect of any sentences in Taylor expansion of a generic function F(z). The 
first term which is a constant will recover a de' Sitter Universe, while the first two terms 
have cover all dynamical properties of the Universe from radiation era to late accelerating 
Universe. No one extra term will affect on the present value of deceleration parameter. 
Adding the third term to the last two terms may make an inflationary behavior, in which 
deceleration parameter may goes back to a huge positive or negative amount. Adding any 
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other terms to the last three terms may change the value of deceleration parameter from an 
unstable value to a stable but negative value in very early Universe. It means that in f(R) 
theory of gravity in the absence of matter, the value of deceleration parameter was not a 
positive value which denotes there was a low inflationary era before main inflation. Adding 
more terms to the last four terms do not change whole the story, but will vary the value of 
deceleration parameter to a larger value in negative direction. Thus for an infinite number 
of terms in Taylor expansion of F(z) the value of deceleration parameter will goes to — oo. 
It means the Universe is born in a main inflationary state. Then if the number of Taylor 
expansion sentences should be finite whether the Universe was born in an inflationary state. 
However the value of expansion coefficients may make some fluctuations between the first 
moment of creation and the beginning of radiation era in f(R) gravity theory. 



[1] Riess, A.G., et al., Observational evidence from supernovae for an accelerating universe and a 
cosmological constant, The Astronomical Journal, vol. 116, no. 3, pp. 1009, 1998. 

[2] Cole, S., et al., The 2dF Galaxy Redshift Survey: Power-spectrum analysis of the final dataset 
and cosmological implications, Mon. Not. Roy. Astron. Soc. 362, 505, 2005 

[3] Spergel, D.N., et al., First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 
Determination of Cosmological Parameters, Astrophys. J. Supple. 148, 175, 2003 

[4] Spergel, D.N., et al., Three- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observa- 
tions: Implications for Cosmology, ApJS 170, 377, 2007 

[5] Nojiri, S., & Odintsov, S.D., Introduction to Modified Gravity and Gravitational Alternative 
for Dark Energy, Int. J. Geom. Meth. Mod. Phys. 4, 115, 2007 

[6] Capozziello, S., & Francaviglia, M., Extended Theories of Gravity and their Cosmological and 
Astrophysical Applications. Gen. Rel. Grav. 40, 357, 2008 

[7] Faraoni, V., f(R) Theories Of Gravity, Rev. Mod. Phys. 82, 451, 2010 

[8] Nojiri, S. and Odintsov, S.D., Phys. Rev. D 74 (2006) 086005, [arXiv: hep-th/0608008], Nojiri, 
S. and Odintsov, S.D., Phys. Lett. B 657 (2007) 238, [arXiv: 0707.1941], Cognola, G., Elizade, 
E., Nojiri, S., Odintsov, S.D., Sebastiani, L. and Zerbini, L., Phys. Rev. D 77 (2008) 046009, 
[arXiv: 0712.4017], Nojiri, S., Odintsov, S.D. and Saez-Gomez, D., [arXiv: 0908.1269], Nojiri, 
S. and Odintsov, S.D., [arXiv: 1011.0544], Elizade, E., Nojiri, S., Odintsov, S.D., Sebastiani, 

9 



L. and Zerbini, L., [arXiv: 1012.2280], 
[9] Amanullah, R., et al., arxiv.org: 1004.1711 



10 



